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Network topology plays an important role in governing the collective dynamics. Partial synchro- 
nization (PaS) on regular networks with a few non-local links is explored. Different PaS patterns 
out of the symmetry breaking are observed for different ways of non-local couplings. The criterion 
for the emergence of PaS is studied. The emergence of PaS is related to the loss of degeneration in 
Lyapunov exponent spectrum. Theoretical and numerical analysis indicate that non-local coupling 
may drastically change the dynamical feature of the network, emphasizing the important topological 
dependence of collective dynamics on complex networks. 
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Synchronization, as a universal cooperative behavior 
and a fundamental mechanism in nature, has been ex- 
tensively studied in relating to numerous phenomena in 
physics, chemistry, and biology 0, S 0- In recent 
years, there has been a growing interest in the synchro- 
nization of spatiotemporal systems, especially in syn- 
chronous dynamics on networks. Synchronization on typ- 
ical complex networks, e^., on small- world networks [5|, 
scale-free networks 0, [Tj or weighted network^ have 
been investigated recently 9|. In spite of these efforts, a 
number of fundamental questions still remain open. For 
example, most of previous works have been focused on 
the onset of global synchronization, while much less was 
explored for partial synchronization prior to the global 
case. However, there exist rich synchronous dynamics 
on complex networks, and many of these behaviors are 
closely related to partially synchronous motions. More- 
over, the mechanism for synchronization on complex net- 
works is still not clear. A good understanding of this 
issue should be relevant to many collective behaviors in 
spatiotemporal systems, especially in complex networks. 

An important feature of complex networks is the exis- 
tence of non-local links. On the one hand, dynamics on 
general complex networks (e.g., small- world or scale- free 
types) are usually difficult to study due to the sophisti- 
cated network topology; On the other hand, simple net- 
works, especially regular networks, cannot be ideal mod- 
els in describing real networks. Therefore, one should 
find some kind of network that is relatively simple but 
complex enough to reflect properties of typical complex 
networks. It is intuitive and feasible to consider a regu- 
lar network plus only a few non-local links, which forms 
the simplest "complex" network. This can be regarded 
as a good bridge between regular networks and general 
complex networks. In this Letter, we address the issue 
of synchronization on rings of coupled oscillators with a 
few shortcut links. We reveal important partial synchro- 
nization (PaS) phenomena induced by symmetry break- 
ing in the presence of non-local couplings. Various syn- 
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FIG. 1: Different Topologies of regular networks with non- 
local links, (a): One shortcut leads to five synchronized 
clusters(IO nodes); (b): One shortcut leads to six clusters 
(10 nodes). (c): One shortcut leads to five clusters (9 nodes), 
(d): Two parallel shortcuts; (e)Two crossed shortcuts; (f)The 
lambda-type network. 

chronous patterns are observed for different types of non- 
local links. This reveals the significant role of network 
topology in governing the global dynamics. We analyze 
in detail the conditions for the emergence of PaS and give 
the criterion of the topological dependence of PaS. 

We apply certain m-dimensional dynamics to nodes in 
a given network: x — f(x),x — (x±, X2, ■ ■ ■ ,x m ). The 
network is formed by using the linear coupling as edges. 
The dynamical equations of the network can be written 
as 

X = F(X) + eT <g> CX, (1) 

where X = (x 1 , x 2 , ■ ■ • , x N ), F = (/\ / 2 , ■ ■ • , f N ), and 
N is the number of nodes, e denotes the coupling 
strength, and T : W n — > R m characterizes the coupling 
scheme among the nodes. C = M — D, where M is the 
adjacency matrix of the network (the element de- 
notes the number of the edges that link node i and j). 
D is a diagonal matrix, and Da = Xw=i-^%- There- 
fore Ylj=i Cij — 0, and its largest eigenvalue Af is (we 
adopt Ai > A 2 > • • • > Xn throughout this Letter). 

We first briefly describe the synchronous phenomena 
we observed in the above networks. Above all, one can 
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define the average distance between two oscillators i and 
j by dij = limy^ 00 (T _1 J \\d% — xj \\dt), and say they are 
in a synchronous state in the sense that dij — > 0. Denote 
the number of non-local links by N s . When N 3 = 1, 
the synchronous patterns are shown in Figs. [Ja)-(c), 
where the nodes with the same symbol will synchronize 
to each other, forming several synchronous clusters. We 
call the formation of these synchronous clusters "partial 
synchronization" (PaS). Several interesting phenomena 
can be observed: (1) PaS emerges not only for oscilla- 
tors directly linked by the non-local couplings, but also 
for those without shortcuts; (2) Synchronous clusters are 
formed by satisfying the mirror symmetry of lattices, i.e., 
if is a pair linked by a shortcut, then (i + k,j — k) 
are all synchronous pairs. More interestingly, if there are 
two non-local links, the emergence of PaS also need not 
only satisfy the mirror symmetry of network, but also 
depend on the linking form. For example, PaS can al- 
ways be achieved on networks with "parallel" links [the 
network with two shortcuts and (i + k,j — k)], as 

shown in Fig. 2Jd). However, for networks with "cross" 
shortcuts [two shortcuts and (i + k,j + k)] shown 

in Fig. nje) and networks with "lambda" -type shortcuts 
[two shortcuts (i,i + k) and (j, j + k)} shown in Fig.^f), 
PaS emerges only for some choices of k. If PaS emerges 
on "cross" and "lambda" networks, synchronization may 
not emerge between pairs that are directly linked. 

The above PaS behavior implies the important role of 
network topology in governing the collective behaviors. It 
is important to understand the mechanism of the above 
interesting PaS phenomena. For example, why does PaS 
appear on networks with mirror symmetry of sites? What 
is the manifestation in the dynamical exponents at the 
onset of PaS? Moreover, can we find a simple criterion in 
exactly determining whether PaS can appear in a given 
network? It is our central task to answer these important 
questions. 

We start with theoretical analysis on the stability of 
complete synchronization |l l| . For the synchronous state 
x 1 ^) = 5? (t) = ■ ■ ■ = x N (t) = s(t), by linearizing Eq. (JTJ 
near the synchronous manifold s(t), i.e., x* 1 = s + u l (i = 
1, 2, • • • , N), one may get U = [I N ®Df+eT®C]U , where 
U = (m 1 , u 2 , ■ ■ ■ , u N ) and Df is the Jacobian function of 
/ near the synchronous state s(t) . Noticing that the first 
term of the new equation is block diagonal with m x m 
blocks, one can diagonalize C. Thus it becomes Vk = 
[Df + e\kT]vk, (k = 1, 2, • • • , AT), where are eigenval- 
ues of C. Then Eq. (Q becomes N independent blocks . 
To judge whether CS can emerge on a network, one needs 
to find the regime where the largest Lyapunov exponent 
of the generic equation v = [Df + aT]v is negative on 
the a plane (a is complex). If all e\k(k — 1, 2, • • • , N) 
are in the negative(stable) regime (except Ai =0), then 
the CS state will be stable. Therefore, the essence of this 
method is to divide the system into two subsystems which 
dominate the dynamics on the synchronous manifold (the 
equation with Ai = 0) and that on the transversal mani- 



fold(the equation with A;(i = 2, 3, • ■ • , N)), and then the 
emergence of the CS is determined by the stability of the 
transversal system. 

For the case of PaS, one also needs the system to be 
divided into two parts, and the stability of one of them 
(the transverse system) determine the emergence of the 
PaS. On the one hand, the symmetry of a network can 
be defined by the invariance under a permutation trans- 
formation. Taking the case in Fig. [Ja) as an example, 
its symmetry requires that F^CF^ 1 = C, where Fn is 
a counteridentity matrix that satisfies i^ ^v—j+i 

= 1, for 

i = 1, 2, • • • , N and F„ = for all i + j ^ N/2 + 1, thus 
the network has a mirror symmetry. It is easy to verify 
that PaS is a solution of the differential equation, then 
each node in a cluster should possess the same dynamics 
|12| . Therefore the mirror symmetry is needed for the 
emergence of PaS. On the other hand, one can introduce 
new variables W = SX, where W = {w 1 , w 2 , ■ • ■ , w N ), 

S = ( \ N/2 F »' 2 ) • (2) 
V In/2 —*n/2 J 

Here In 12 IS a ^/2-order identity matrix. For a network 
with a mirror symmetry, the above transformation can 
divide the system into two subsystems. Although this 
requires the adjacent matrix C to be persymmetric (i.e., 
C is symmetric about its anti-diagonal), the condition 
can always be satisfied by relabelling the node for a net- 
work with a mirror symmetry. Furthermore, if a network 
satisfies the mirror symmetry and its label obeys the rule 
in Fig. QJ, the adjacent matrix C can be block diagonal- 
ized by the transformation J2J: 

This transformation only re-distributes the eigenvalue 
spectrum of the matrix C into two blocks A and B. Since 

Sj=i A%j = 0) the largest eigenvalue of A is 0. 

By using the above transformation and considering 
that we only care about the dynamics near the syn- 
chronous manifold, the system can be separated into two 
parts, and they can be linearized near their synchronous 
manifold respectively: 

SW im = [DF{W m ) + eTA]SW im , (3) 

5W tv = [DF(W tv ) + eTB]SW tv , (4) 

where W im = (w 1 , w 2 , ■ ■ ■ , w N / 2 ) and W tv = 
(uj N / 2+1 ,iv N / 2 + 2 • • ,w N ). For the PaS case discussed 
above, we have Wi m = X\ + X 2 ,Wtv = X\ — 
X 2 , where Xi = (x 1 , x 2 , ■ ■ ■ , x N l 2 ), and X 2 = 
{x N / 2+l ,x N / 2+2 ,--- ,x N ). Here, W im denotes the dy- 
namics on synchronous manifold and Wtv denotes the 
dynamics on transversal manifold. 

Back to Eqs.© and J3J. If one increase the coupling 
from 0, synchronization will be achieved for (3) and (4) 
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FIG. 2: (a): The label rule for N s = 2; (b)-(d): The compar- 
ison between e c - e Pa s and XT - X\ v , N=20. XT - X\ v are 
amplified by a proper times; (b): Parallel network, j=3; (c): 
Cross network, k=9; (d): Lambda network, j=3; (e): Topol- 
ogy graph of the network for N s = 3; (f): The same as (b)-(d) 
for N s = 3. 

when eA™ (AJ 71 is the second largest eigenvalue of A) and 
eX\ v (X\ v is the largest eigenvalue of B) respectively fall 
in the stable regime of the a plane. The synchronization 
(the stability of synchronous state) of Eq. J3J suggest the 
emergence of PaS. If A™ < Aj"", the synchronization of 
Eq.© will be achieved at a smaller coupling than that 
of Eq.lQJ, thus a global synchronization state has been 
achieved when Eq.(0} is stable. Therefore the emergence 
of PaS requires 

A| m > Af . (5) 

Now let us test the validity of this criterion and numer- 
ically explore the PaS phenomena. We take the Lorenz 
oscillator as the node of the network: dx/dt = o(y — x), 
dy/dt — rx — y — xz, dz/dt — xy — bz. Here a = 10, r = 
27, b = 8/3. We adopt the coupling scheme by the 3x3 
matrix T as P^i = 1 and Tij = for i ^ 2,j ^ 1. Let 
us first try to find all possible networks that satisfy the 
mirror symmetry. Taking the network with two non-local 
links for instance (the case of N c = 1 can always topo- 
logically lead to PaS). For convenience we always take 
node i = 1 as one of the ends. The second and the third 
ends are labelled as j and k (k > j), respectively. Due 
to the mirror symmetry, the label of the fourth end is 
714 = k + j — 1. There are three possible links that hold 
the mirror symmetry , as shown in Figs.^d),(e) and (f). 
Noticing the translational symmetry of the ring network, 
we have 2 < j < N/2 and j < k < N/2. Thus by vary- 
ing j and k, all symmetric networks with two non-local 
couplings can be found [Fig.[2{a)]. 

We use Spas and ec to represent the critical coupling 
of PaS and CS, respectively. Figs^b) and (c) give a 
comparison between differences AA = X 1 ™ — X^ and 
Ae = ec — epaS f° r different networks. It is instruc- 
tive to find that AA passes zero when Ae becomes zero. 
Ae — > means the coincidence between CS and PaS, 



implying that PaS doesn't appear. This good correspon- 
dence verifies our criterion in JSJ). Our investigations re- 
veal that one-shortcut networks and "parallel" networks 
can always achieve PaS because both AA and Ae are 
positive for all k, as shown in Fig[3fb) for the "parallel" 
case. For "cross" and "lambda" networks, PaS occurs 
only for some configurations, as shown in FigsEJc) and 
(d). To illustrate these results, we plot the PaS critical 
line AA = in the (2j/N,2k/N) plane. FigEIa) plots 
the PaS regime for "cross" networks, and Fig|^b) gives 
the phase diagram of "lambda" networks. Different lines 
correspond to different nodes number. These lines co- 
incide with each other pretty well, indicating that PaS 
depends on the topology of networks and is independent 
of numbers of nodes. 

The above criterion for PaS can also be applied to net- 
works with more than two non-local links, while the dis- 
cussions become more complicated. However, the above 
phenomena and criterion still remain valid. Fig|5jf) 
shows one case of three shortcuts. There are 20 nodes 
in a network, and two non-local links are fixed to be 
(11,3) and (11,19), and the third link is (ll+m,ll-m), 
where m = 1, 2, • • ■ ,9, see Fig|!Je). PaS occurs only for 
the m = 4, 5 cases, as shown in Fig^f). 

Let's observe the manifestation of the Lyapunov expo- 
nent (LE) spectrum (LES) of the network at the onset of 
PaS. FigOJa) gives the first four largest LE's for a ring of 
6 nodes without shortcuts. One can find that the third 
and the fourth largest LE's pass zero at the same cou- 
pling (e = ec), indicating the degeneration of LE's at the 
onset of CS for ring network. Fig^b) is the four largest 
LE's of a ring network with one non-local link (1,4). Ob- 
viously the presence of non-local links leads to a break- 
ing of degeneration of LES. The fourth largest LE passes 
zero at a very small coupling. It also can be found that 
^1,41^2,3)^5,6 become at the same coupling, as shown 
in Fig0Jc), but ^1,2,^3,4)^1,6)^4,5 / 0. Therefore, PaS 
emerges if and only if there only are two positive LE's. 
This conclusion was tested to be true for different node 
numbers and different shortcut numbers and link ways. 

These results can be interpreted by analyzing eigen- 
values of the adjacent matrix. The adjacent matrix of 
a ring network is a circulant matrix, its eigenvalues are 
conjugated besides the largest and the smallest ones. So 
when the coupling increases, eA2 and eA3 will fall into 
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FIG. 4: (a): The first four largest LEs of a ring network; (b): 
The first four largest LE's of a ring network with one shortcut 
(1,4); (c): average trajectory distances di,4, (£2,3, ds,6 (solid 
fines), and rfi,2, ^3,4, di,e, ^4,5 (dashed lines) varying with the 
coupling, where the network is the same as in (b); (d): The 
three largest LE's of the synchronous manifold (solid lines) 
and the largest LE of the transversal manifold. 

the stable regime at same critical coupling, above which 
the third and fourth LE's corresponding to A2 and A3 be- 
come negative (there are always at least one positive LE 
and one zero LE for a chaotic system). When non-local 
couplings are applied, the degeneration in the eigenval- 
ues and LES is lost. Since the emergence of PaS requires 
A| m > A* 11 , when eX^" is in stable regime, there are only 
two eigenvalues (A™ = and A?> m ) outside of the sta- 
ble regime. That is, there are only two positive LE's. 
One can also calculate the conditional LE spectrum of 
the synchronous manifold (Eq.lJSJ) and its correspond- 
ing transversal manifold (Eq.Q). It can be found from 



Fig^Jd) that there are only two positive conditional LE's 
of synchronous state(solid line) when the largest condi- 
tional LE of the transversal manifold (dotted line) passes 
zero. 

In conclusion, in this letter we explored synchroniza- 
tion behaviors of complex networks with a few non-local 
links. The criterion for the emergence of PaS is the- 
oretically given and verified in numerical simulations. 
The manifestation of Lyapunov exponent spectrum is 
discussed for PaS. Our results and criterion proposed in 
this Letter are general, and they are independent of the 
dynamics on sites and coupling forms, we tested differ- 
ent dynamics on the nodes, such as Lorenz oscillators, 
Rossler systems, and nonlinear maps (e.g., the logistic 
map and circle maps). All these numerical experiments 
support our criterions and conclusions. These "simple" 
complex networks are not only simple but also typical, 
and it can be a useful tool to study more complex net- 
works. The synchronization behavior we studied in this 
letter reveal the importance of network topology. Differ- 
ent types of shortcuts may drastically change the topol- 
ogy and furthermore the dynamics of a network. Our 
work should illuminate a good step in understanding the 
synchronization processes on general complex networks 
from a dynamical point of view. 
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